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oo- 1 Introduction 

In |NV09j . I. Nourdin and F.G. Viens have introduced sufficient conditions to prove the existence 
of a density for a Malliavin differentiable random variable and to provide upper and lower Gaussian 
estimates for this density. 

' This result has led to several research papers, such as those by D. Nualart and L. Quer-Sardanyons 

■ ( |NQ09a| , |NQ09b| ) , in which these authors applied Nourdin and Viens result to solutions of quasi- 

linear stochastic partial differential equations and to a class of stochastic equations with additive 
noise. 

PsJ ■ In this paper, we use Nourdin and Viens's approach to prove that, under proper conditions on 

<*' . the coefficients, each component of the solution (Xt, Yt, Zt) to a backward stochastic differential 

QQ I equation 

m 
in 

o: 



Xt = xo+ [ h{Xs)ds+ [ a{Xs)dWs (1.1) 
Jo Jo 

Yt = i+ f f{Xs,Ys,Zs)ds- [ ZsdWs (1.2) 
Jt Jt 



has a density for which upper and lower Gaussian bounds can be derived. This implies to study 
r> . the relation between the coefficients of the diffusion equation (|1.1|) and the coefficients of the 

C^- backward SDE ([12]). 



Our paper is organized as follows: in a first part, we study the component Yt of the solution and 
in a second and last part, we focus on the component Zt (for which, up to our knowledge, no 
density existence results exist). We will not develop a specific study of the first component Xt 
of the solution to the BSDE due to the fact that the question of the existence of a density for 
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the solution to an SDE of the type (jl.ip and the properties of this density has been intensively 
studied and we refer the reader to |Nua06j for an extensive survey of the existing litterature and 
results on this topic. 

Equations of the type (|1.2p were introduced in |PP90j and are closely related with viscosity 
solution to PDEs. These equations have been intensively studied and have many applications in 
control theory and financial methods among others. 

The existence of the density for the random variable It at a fixed time t € (0, T), as well as upper 
bounds for its tail behavior, have been proven by F. Antonelli and A. Kohatsu-Higa |AK05j . using 
the Bouleau-Hirsch Theorem. We retrieve Antonelli and Kohatsu-Higa's existence result for the 
density of Yt, and we also derive Gaussian estimates for it. In order to provide (additionally to 
the existence result itself) estimates for the density of 1^, we need to strengthen the hypotheses 
of Antonelli and Kohatsu-Higa. 

We also address the question of the existence of a density for the random variable Zt as well as 
the possibility of deriving Gaussian estimates for it. This question has not been solved in |AK05j . 
We need the same hypotheses as in the case of Yt, as well as additional ones, since Zt can be 
expressed as a function of the Malliavin derivative of If. 

In order to be self contained, we first give an overview of some elements of Malliavin calculus in 
Section 2, and the corresponding notations. Section 3 is dedicated to the component Yt of the 
solution to the BSDE. Section 4 deals with the question of the existence of a density for Zt and is 
organized in two subsections, dealing respectively with the question of the existence of a density 
and conditions for this density to be bounded by Gaussian upper and lower estimates. 

2 Framework, main tools and notations 
2.1 Elements of Malliavin calculus 

Consider the real separable Hilbert space L^{[0,T]) and let (W{ip),ip£ L'^{[0,T])) be an isonormal 
Gaussian process on a probability space (0, 21, P), that is, a centered Gaussian family of random 
variables such that E {W {ip)W {^p)) = {f,tp)L^([o,T])- For any integer n > 1, denote by /„ the 
multiple stochastic integral with respect to W (see |Nua06j for an extensive survey on Malliavin 
calculus). The map is actually an isometry between the Hilbert space L^([0,T]") equipped 
with the scaled norm -^\\ ■ \\l^{[o,t]") and the Wiener chaos of order n, which is defined as the 
closed linear span of the random variables Hn{W{ip)) where if € L^([0,T]), ||<^||l2([q = 1 and 
Hn is the Hermite polynomial of degree n > 1, that is defined by 




The isometry of multiple integrals can be written as follows: for positive integers m, n 



'^K/>f)L2{[o,T]") ifm = n, 
if m 7^ n. 



It also holds that 



In{f)=In{f) 
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where / denotes the symmetrization of / defined by 

f{xi,...,Xn) = X] /(^fT(l), • • • ,2;f7(n))- 

We recall that any square integrable random variable F which is measurable with respect to 
the a-algebra generated by W can be expanded into an orthogonal sum of multiple stochastic 
integrals 

F = J2Ufn) (2.1) 

n>0 

where /„ G L^([0,T]"') are (uniquely determined) symmetric functions and /o(/o) = E [F]. 

Let L be the Ornstein-Uhlenbeck operator defined by LF = — X]„>o "-^n(/n) if F is given by 
(j2.ip and satisfies X^„>i ?T-^?T-!||/n|P < oo. For p > 1 and a G M we introduce the Sobolev- 
Watanabe space B"'^ as the closure of the set of random variables of the form ()2.2p (see (1-28) 
in [Nua06j ) with respect to the norm defined by 

||F|U,p = ||(/-L)tF||ip(a), 

where / represents the identity. We denote by D the Malliavin derivative operator that acts on 
smooth random variables of the form 

F = g{W{vi),...,W{ipn)), (2.2) 
where 5 is a smooth function with compact support and ipi G L^([0,T]), as follows: 

n ^ 

DF = Y,Q^^{W{^^),...,W{ipn))^i. 

i=l 

The operator D is continuous from D"'^ into B"~^'*' (L^([0,T])) . The adjoint of D is denoted 
by S and is called the divergence (or Skorohod) integral. It is a continuous operator from 
pa,p ^^2qq^ T])^ into 0"+-*^'^. More generally, we can introduce iterated weak derivatives of order 
k. If F is a smooth random variables and A; is a positive integer, we set 

A',...,t,^ = AiA....A,F. 
We have the following duality relationship between D and 6 for F G D^'^ and u G dom 6 

E(F5(n)) = E(Z)F,n)i2([o,T]). 

For adapted integrands, the divergence integral coincides with the classical Ito integral. We will 
use the notation 

6{u) = / UsdWs- 
Jo 

Note that the following integration by parts relation between D and 6 holds 

Dt{6{u)) =ut+ [ DtUsdWs, 
Jo 

where u G B^'2(l2([0, T])) is such that 6{u) G B^'^. 
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2.2 Density existence and Gaussian estimates 

A classical density existence result is the celebrated Bouleau-Hirsch theorem (see |Nua06j for an 
extensive survey on this result). This result provides conditions in term of Malliavin derivatives 
for a random variable to have a density. 

Theorem 2.1 (Bouleau-Hirsch). Let F be a random variable of the space D"^'^ and suppose that 
ll-^-^llL2([o,r]) > Then the law of F is absolutely continuous with respect to the Lebesgue 
measure on M. 

In |NV09j ■ Corollary 3.5, Nourdin and Viens have given the following sufficient condition for a 
weakly differentiable random variable to have a density with lower and upper Gaussian estimates. 

Proposition 2.2. Let F be in B^'^ and let the function g be defined for a// x G M by 

g{x) = E [{DF, -Z)L-iF)i2([o,T]) " E(F) = x) . (2.3) 
// there exist positive constants 7mirn7max such that, for all x G M, almost surely 

< 7min ^ 9{x) ^ 7max 

then F has a density p satisfying, for almost all z & R 

E|F-E(F)| / (z-E(F))2\ E\F-E{F)\ ( [z-E[F)f 



72 ^2 17-2 exp. — 



2 7^ V 2y . / - / V / - 2 y . ' V 2^' 

/max \ ^ /mm / ^ 'mm V 'max 



Furthermore, Nourdin and Viens have also provided the following useful result, which gives some 
rather explicit description of g{x). Recall that W = [W{(j)),(p £ ([0, T])). 

Proposition 2.3. Let F be in 9^'"^ and write DF = ^p{W) with a measurable function : 
]^l2([o,t]) _^ L2([0,r]). Then, if g{x) is defined by (EJ), we have 

g{x) = ^"e""E (E'((<I.^(t^),$|(VF))i2([o,T]))|i" - E(F) = x) du, 

where ^'^{W) = ^F{e~^W + Vl — e~'^^W'), W' stands for an independent copy of W , and is 
such that W and W are defined on the product probability space (O x fi', T' ,¥ x ¥') and E' 
denotes the mathematical expectation with respect to P'. 

2.3 Notations 

We denote by C^(M^) the space of n-times differentiable functions on W with bounded partial 
derivatives up to order n. 



Let / be a three times differentiable function of three variables x, y and z. We will use the 
following notations : §| = = /y, || = g = /^^, ^{x,y) = fyy, £^ = f^y, 



fyx- 



dydx J ' 

We will also use the following notation for the Lie bracket : [h, g] = hg' —gh' , where h, g 
In the whole paper, c and C will denote constants that may vary from line to line. 
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3 Density of Yf : existence and Gaussian estimates 



The following backward stochastic differential equation was introduced in Pardoux and Peng 
|PP90j (see also |PP92j ) and density properties of its solutions were investigated in |AK05] : 



Xt = xo + 



f b{Xs)ds+ [ a{Xs)dW, 
Jo 

Yt = C+ f fiXs,Ys,Zs)ds- [ Zs 
Jt Jt 



(3.1) 
(3.2) 



In this section, we give conditions for the random variable Yt to have a density which can be 
bounded from above and below by Gaussian ones. 



3.1 Hypotheses 

We consider 6, a and / to be appropriately smooth functions to ensure the existence and unique- 
ness of solutions to equations ()3.ip and ()3.2p . We also impose additional conditions needed to 
state our main result: 

' HI : e G i^(^^,-7^T)nB^'2 and V < T, < c < < C a.s. 
< H2 : / G Cfe^(M^) and < < C. 

^ H3 : 6 G Cl{M.), a £ C^iR), 0<a<C and < Ma. 

where [b, a] denotes the Lie bracket between b and a. 

Remark 3.1. Note that it is natural to have a condition on the Lie bracket between b and a 
as this quantity is the one that appearing in the classical density and smoothness results (e.g. 
Hormander's brackets condition). 

Remark 3.2. The hypotheses on the positivity of a and fx are made in order to make the proofs 
as readable as possible. In fact, one only needs a and fx to have the same sign to draw to same 
conclusions. 



3.2 Main result 

Under the above assumptions, Yt has a density for which the following Gaussian estimates can 
be derived. 

Theorem 3.3. Under the above hypotheses, for t € (0,T) the random variable Yt defined in 
113. g|) has a density pYt ■ Furthermore, there exist strictly positive constants c and C such that, for 
almost all y & M and all t € [0,T], py^ satisfies the following: 

E\Yt-E{Yt)\ ( {y-nYt)? \ 

^ p'-^^y^ ^ — 2ci — [ — 2^t — J • 

The rest of this section is devoted to the proof of Theorem 13.31 which is divided in three steps. 
In the first step, we will prove that the Malliavin derivative of Xt is bounded and non-negative. 
This property will be necessary in Step 2 (as DXt appears in the Malliavin derivative of Yt). 



E\Yt-E{ Yt)\ 
2ct 



-exp 



(j/-E(r,))- 
2Ct 
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In the second step, we will derive upper and lower bounds for the Malliavin derivative of Yt. 
Indeed, our purpose is to use the Nourdin-Viens formula (Proposition I2.2p in which one needs 
to bound a function of the Malliavin derivative of random variable for which density results are 
investigated. We will make use of the properties of DXt prooved in Step 1. 

In the third step, we will use the bounds obtained on DYt and the Nourdin-Viens formula to 
conclude the proof. 



Proof of Theorem [3T3t 



Step 1 : Boundedness and positivity of DXt 



Consider equation (j3.ip . Using a Lamperti transform (see |Lam64j or |KS91j pp. 294-295 exer- 
cise 2.20), we compute the Malliavin derivative of Xt. The Lamperti transform of Xt, hereafter 
denoted by Ut, is given by 



where 



Ut = g{xo)+ [' Pog~\Us)ds + Wt, 
Jo 



9{x) = I and /3(x) = -(x) - 

0'(iij cr 2 



Computing the Malliavin derivative of Ut yields, for 9 G [0,t], 

DeUt = 1 + A/3 o g-^)'{Us)DeUsds = exp [ / 
Je VJe 



{Pog-^)'{Us)ds 



(3.3) 



Deriving the identity gog ^(x) = x on g{supp{Xt)) yields (g ^)'(x) = ao g ^(x). Using this fact 
we get (/3 o g~^y{x) = /3' o g~^{x){g~^y{x) = (/3'o") o g~^[x). In addition, it is easy to check that 
on g{supp{Xt)), 



(/3V)(x) 



[a,b]{x) {aa"){x) 



cr(x) 



(3.4) 



Gathering those results and using hypothesis (H3) of Subsection 13.11 immediately yields on 
g{supp{Xt)) 

-C < (/3oG-i)' < C7, 
where C is a positive constant. Using (j3.3p . we deduce, P-a.s, 

< c < DgUt < C. (3.5) 
Furthermore, as Xt = g~^{Ut), it holds that, for < < t < T, 

DeXt = {g-^y{Ut)DeUt = a o g'\Ut)DeUt. (3.6) 
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Combining (j3.5p and ()3.6p with the fact that a is bounded and non-negative yields, P-a.s, 

< DeXt < C. (3.7) 



Step 2: Computation of bounds on DYt 



We at first represent -Dglt by means of an equivalent probability; this is similar to |AK05j and 
the proof is included for the sake of completeness. It is well known (see for example Theorem 2.2 
in |AK05j ^ that, for every t £ (0,r], Yt G D^'^ and Z G {0,T;O^'^). Furthermore, since 9 < t, 
we have 

DeYt=De^-j^ DeZ^dWs 

+ [U (X,, Ys, Zs) DeXs + fy (Xs, n, Zs) DgY, + /, (X„ y„ Z,) DgZ,] ds. (3.8) 

The product e^o fy'^^-^-^^-)'^' DeYt yields a more suitable representation of DgYt; indeed, for t G 
(0,r], and < 6* < t 

d Uo fy{Xs,Ys,Zs)ds-^^Yt] = iDeYte^O fy(Xs,Ys,Zs)ds ^^^^ Zt) 

_e/o fy{x.,Ys,Zs)ds (^^^ Yt, Zt) DeXt + fy {Xt, Yt, Zt) DeYt 
+/, {Xt,Yt, Zt) DeZt)] dt + e/o fy^^'^Y^^^^)^^ DeZtdWf 
Integrating from t to T yields, for 9 < t, 

/T 
e/o fy{Xr,Yr,Zr)dr (X„ y„ Zs) DeXs 

j-T 

+/, {Xs,Ys,Zs) DgZs] ds + e/o M^'-'^'-'^^^'^- DeZ^dWs. 
Note that DqYt = DqC; therefore, for t G (0,T], 

Del^t =e/,^ fy(x.,Ys,z.)dsj^^^ ^ /,(X.,y,.,Z.)dr (^^^ ^^^^ 

/■T 

+/, (X„ n, Z,) - e-lt fy^^^'^^^^^^'^'DeZsdWs. 

Let Wt = Wt — /q (Xs, Yg, Zs) ds. Because fz G C^(M), Novikov's condition is verified and W 
is a Brownian motion under some equivalent probability P. Girsanov's theorem yields 

DeYt =e^^ fyiXs,Y.,Zs)dsj^^^ ^ ^J/ /«(x,,y,,z.)dr^^ (X„ n, Z,) i^eX.ds 
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Conditionning by Tt under P and taking into account the fact that and DqYi are adapted with 
respect to Tt while fy (X^, 1^, Zr) dr and DqZs are J-"s-adapted for 6 < t < s <T, we obtain 



DgYt =E ( e-Zi fy(Xs,Ys,Zs)dsjj^^ 



+ E 



(3.9) 



Using hypotheses (HI) and (H2), the first summand in ()3.9p can be bounded by two positive 
constants c and C in the fohowing manner: 

'""^ TA<C. (3.10) 



Using the results on DXt proved in Step 1 along with hypothesis (HI), we deduce that the second 
summand in ()3.9p is bounded and non-negative: there is a positive constant C such that 

< E fyi^r.Yr,Zr)dr f^^x^^Ys,Zs)DeXsds < C. (3.11) 

Combining the bounds p.lOp and ()3.1ip . we immediately deduce that there exist two positive 
constants c and C such that 



< c < DeYt < C. 



(3.12) 



Step 3: Conclusion of the proof by the Nourdin Viens formula 

Write D,Yt = ^y^{W) with a measurable function : M^'([0'^]) — y L2([0,T]). Then the 
bounds obtained in (j3.12p yield, for 9 < t, 

< c < $y,(VF) < C. 



Define ^y^iW) = $^^(e-"l^ + VT^^e^W) for u G [0,+oo[, where W stands for an inde- 
pendent copy of W and is such that W and W are defined on the product probability space 
(0 X T' , P X P'). It is clear that, for 6 <t,we have for any u G [0, oo), 

< c < ^^Y^{W) < C. 
Combining the two previous bounds yields, for ^ < t, n G [0, oo), 

Q<c^ <¥y^{W)^^^^{W)<C'^. (3.13) 
Using the notation from Propositions 12.21 and 12.31 we define 

9{y) = e-«E l^j\'y^{W)^'^';{W)de^ \Yt - ^{Yt) = du. 

The bounds obtained in (|3.13p immediately yield 

< ct < g{y) < Ct, 

with strictly positive constants c and C. Thus, Propositions 12.21 and 12.31 conclude the proof of 
Theorem [331 □ 
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4 Density of Zt : existence and Gaussian estimates 



We consider equations (j3.1|) and (|3.2|) with a function /* that only has a Hnear dependency on 
Z, i.e. 



where z) = /(x, y) + qz, a G M. 

Because of the dependency of / on Z, the Malhavin derivative DZ will depend on D^Z^ which 
is not suitable for analyzing it within our framework. One can circumvent the above mentioned 
issue by using the Girsanov theorem to dispose of the impeding terms (similarly as done in the 
proof of Theorem 13. 3|) . In order to clarify the proofs and to improve readability, we will consider 
that this step has already been performed in all of our proofs. This procedure leaves us with an 
equation of the type 



which is the one that will be referred to in the proofs of the upcoming results. 

In the following subsections, we will prove our two main results concerning Zt- We will be- 
gin by giving sufficient conditions for Zt to have a density. Up to our knowledge, this is the first 
result on density existence for the component Z of the solution to equation ()4.2p . We will then 
study in what framework and under what conditions this density can be bounded by Gaussian 
estimates. 

4.1 Existence of a density for Zt 

We list in the next section the full set of hypotheses we need in this section. 
4.1.1 Hypotheses 

We consider 6, a and /* to be appropriately smooth functions to ensure the existence and unique- 
ness of solutions to equations ()4.ip and ()4.2p . We also impose additional conditions needed to 
prove Theorem 14.31 




(4.1) 



(4.2) 




(4.3) 



(4.4) 



' H4 : ^ G L^{n,TT) D D^'^, W < T, Dei > a.s and V 6* < t < T, Z?| i^ > a.s. 
< H5 : / G C^(M2) and fx, fy, fxy, fxx, fyy > a.s. 

, H6 : be C2(M), a G C^(M), a, a', -a", -a'" > a.s and [cj, [a, h]] > a.s. 



where [b, a] denotes the Lie bracket between b and a. 



9 



Remark 4.1. Note that it is natural to have a condition on the iterated Lie bracket [o", [c, 6]] 
between b and a as second order MaUiavin derivatives appear in the expression of Zt. 

Remark 4.2. The hypotheses on the signs of a and fx are made in order to make the proofs as 
readable as possible. It is possible to have more complex hypotheses for the signs of the products 
of derivatives of a and derivatives of /. 

The following theorem states that under the above hypotheses, Zt has a density on M 

Theorem 4.3. Under the above hypotheses, for t G (0, T) the law of the random variable Zt is 
absolutely continuous with respect to the Lebesgue measure on M. 

Before proving Theorem 14.31 we will first give a technical Lemma and a Proposition which will 
play a key role in the upcoming proof of this Theorem. First recall a lemma used to calculate 
the Malliavin derivative of a product of random variables in D^'^ (for example, see |Nua06j . p. 36, 
exercice 1.2.12). 

Lemma 4.4. (i) Let s, t £ [0, T] and F G B^'^ . ^/^g^ /^^^g ^ {F\Tt) G O^'^ and 

DsE{F\Tt) = 'E^{D,F\Tt)ls<t. 
(a) If F,G £ B^'^ are such that F and ||^^||l2([q 'pj) are bounded, then FG G B^'^ and 

D{FG) = FDG + GDF. 

The rest of this section is devoted to the proof of Theorem 14.31 which is divided in three steps. 
In the first step, we will prove that under the conditions of subsection 14.1.11 the second-order 
Malliavin derivatives of X and Y are non-negative. This will be of importance in Step 2 (as these 
second-order derivatives appear in the expression of DZt- 

In the second step, we will show that DZt is positive almost surely. This will ensure that 
\\DZt\\L-i{[o,T\) > a.s. 

In the third and last step, we will use the Bouleau-Hirsch Theorem to conclude the proof (see 
Theorem 12. ip . 

Proof of Theorem 14. 3t 

Step 1: Non-negativity of D^X and D^Y 

We start by proving that for 0<^<t<s<r, P-a.s Dg ^Xg is non-negative. 

Applying the Malliavin derivative to (|3.6p and using the second point in Lemma 14. 4t we deduce 
for 6*, i < s < T, since Us = g{Xs), 

DltXs =(f7 o g~^)'{Us)DeUsDtUs + {o o g-^){Us)DltUs 

={a'a){Xs)DeUsDtUs + a{Xs)DltUs. (4.5) 

Hypothesis (H6) ensures that the term [a' a){X s)DqU sDtU s is non-negative. It remains to prove 
that the second summand in (|4.5p is also non-negative. As a is non-negative, we focus on proving 
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that Dq^Us is too. Applying once again the Mahiavin derivative operator to (|3.3p and using the 
second point in Lemma |4.4| we deduce for 9 < t < s, 

Dl ^Us = j\p o g~y{Ur)DtUrDeUrdr + J\(3 o g'^)' {Ur)DltUrdr 

= J' e^riP°9-')'{u.)dv^p ^ g-^)" {Ur)DtUrDeUrdr 

(/3 o g-^)"{Ur)DrUsDtUrDeUrdr. 
Further calculations yield the following expression 

[f3og } (x) = lal— -^j - -{a a) aj og {x) 

f[a, [a,b]] 1 / „ N/ \ X 
= ( ^{^ V '^]°9 (x). 

Again, Hypothesis (H6) ensures that the term (/3 o g~^)" (x) is non-negative and thus that Dg ^Ug 
is too. We have finally obtained that D'^JJs > a.s. 

We will now focus on Dq ^Yg and prove that for 0<9<t<s<T, P-a.s it is also non- 
negative. Applying once more the Malliavin derivative operator to DqYs in ()3.8p and using the 
second point in Lemma 14.41 since / does not depend on Z we obtain, for 0<9<t<s<T, 

Dl,Ys=Dl,C- r Dl.ZrdWr 

J s 

+ {Xr, Yr) DeXrDtXr + U i^r, Yr) Dj ^Xr 

+ fy., {Xr, Yr) {DeYrDtXr + DgXrDtYr) 
+ fyy {Xr,Yr) DeYrDtYr + fy {Xr, Yr) OltYr^dr. 

Since Dg ^Ys solves a linear equation and is J-'^-measurable, we have that, for < ^ < t < s < T, 
Dl,Ys=B[eI^M^^'^^^'^Dj4\T?l 

+ MXu,Yu)du^j^^ (^^^ Yr) DeXrDtXr + {Xr, Yr) Dj.Xr 

+ fy^ {Xr,Yr) {DeYrDtXr + DeXrDtYr) + fyy {Xr, Yr) DeYrDtYr^dr^Ts) . (4.6) 

Using hypotheses (H4) and (H5) along with the fact that forO<9<t<s<T, P-a.s Dj^Xs 
is non-negative, we obtain that for 0<9<t<s<T, P-a.s Dg ^Y^ is non-negative. 

Step 2: Positivity of DZt 

Using a representation of Z, we compute its Malliavin derivative and show that under the hy- 
potheses of Subsection 14.1.11 it is almost surely positive. We begin by giving a representation of 
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Z. We do not use the one from jPP92j in terms of gradient, that is Zt = cf {Xt) (VX^)^^ VYj, 
but rather use the fact that Zt can be represented by use of the Clark-Ocone formula. Indeed, 
by the uniqueness of the solution (y, Z), Zt can be written as 



Zt = B[Dt^ + Dt £ fiXs,Ys)ds ^ ^ -^1-2 



J^t G E)^'^. (4.7) 



Using this fact, we get for t e [0,T] 

Zt =E (^Ae + {fx {Xs,Ys) DtXs + fy {Xs,Ys) DtYs} ds\Tt^ . 

Let 9 < t. We use both points of Lemma 14.41 in order to calculate the first order Malliavin 
derivative of Zt. This leads, for 6 < t: 

DeZt =^(pl^i + {Xs.Ys) DeXsDtXs + fy,, (X„ Y^) {DeYsDtXs + DeX^DtYs) 



+ fyy (Xs, Ys) DeYsDtYs + f; (X„ y,) Dj ^Xs + fy (X„ y,) Dl^Ys]ds 

Using Hypotheses (H4) and (H5) along with the results obtained in Step 1, we obtain that for 
0<6'<t<r, P - a.s, DgZt > 0. 

Step 3: Conclusion of the proof by the Bouleau Hirsch Theorem 

For all t < T, we have 

T 



l-^^t|lL2([0,T]) = {DgZtfde 



Using the fact that for < < t < T, P — a.s, DgZt > proved in Step 2, we deduce that 
||DZi||£,2QQ > a.s. Applying the Bouleau-Hirsch Theorem (see Theorem 12. ip concludes the 
proof. □ 

Remark 4.5. Theorem 14.31 has been proven under a set of hypotheses (those of Subsection 14. 1 . 1 1 ) 
based on the fact that a is positive. The case where a is negative was included neither in the 
proof nor in the hypotheses for the sake of clarity and readability of the paper. However, as 
already mentioned in Remark 14. 2 1 this case can be addressed (without any further difficulties) by 
using the following transformations: o" — )• fi := —a and W ^ W := —W. After performing those 
tranformations, it suffices to consider {X , Y, Z) = {X, Y, —Z) to be the solution of 



Xt = xo+ [ b{Xs)ds+ [ a{Xs)dWs 
Jo Jo 

Yt = C+ [ f(Xr,Yr)dr- [ ZrdWr 



This brings the problem back to the set of hypotheses of Subsection 14.1.11 and it can be dealt 
with using the techniques presented in the last section. 
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4.2 Gaussian bounds for the density of Zt 

In this section, we study a particular case of equations ()4.3p and ()4.4p and show that under proper 
conditions, the density of Zt can be bounded from above and below by Gaussian estimates. The 
backward equation we study is the following: 

= {Wt.) + /(y,)ds - Z,dW,, (4.9) 

4.2.1 Hypotheses 

We consider / to be an appropriately smooth function to ensure the existence and uniqueness of 
solutions to equation ()4.9p . We also impose additional conditions needed to prove Theorem 14.61 

H7 : (/. G C^(M) and 0" > c> 0. 
H8 : / G C2(M) and /,/">0. 

The following theorem states that under the above hypotheses, Zt has a density that can be 
bounded from above and below by Gaussian estimates. 

Theorem 4.6. Under the above hypotheses, for t G (0, T) the random variable Zt defined in 
i4-^ a density pzt ■ Furthermore, there exist strictly positive constants c and C such that, 
for almost all y G M, pz^ satisfies the following: 

E|Zt-E(Zt)| / (z-E(Zt))2\ , , E|Zt-E(Zt)| [ {z-^{Zt)f 
-exp iTTT, < PzAz) < ^t;: exp 



2ct ' \ 2Ct J - r^ty , - 2Ct V 2ct 

Proof: We will proceed in two steps, the first one being dedicated to proving that the Malliavin 
derivative of Zt is bounded and bigger than a positive constant. The second step will be to use 
the Nourdin-Viens formula to conclude the proof. 

Step 1 : Boundedness and positivity of DZt 

Yt being defined as in equation (|4.9p . equation ()3.9p becomes 

DeYt = E (ei'*^^'(^^)'^>'(iyT)| -^t) , 

and equation ()4.6p becomes 

Dl,Y, =B(^ef^ f'(^^^''^<l>"{WT)\Ts) +B(^l\lsf'(y^)duf"(Y,)DeYrDtYrdr\Ts). 

Using hypotheses (H7) and (H8), we obtain that < \DgYt\ < C and < c < Dj^Ys < C. We 
finally compute DgZt from equation ()4.8p and we get 

DeZt =b(4>"{Wt) + [f"{Ys)DeY,DtYs + f (Y,) dI^Y^} ds\Tt) . 

Using Hypotheses (H7) and (H8) again, we finally get 

< c < DgZt < C. (4.10) 
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Step 2: Conclusion of the proof by the Nourdin Viens formula 

Write D,Zt = ^z^{W) with a measurable function : M^'([0'^1) — > L^{[0,T]). Then 
yields, for < t, < c < <^%{W) < C. As previously done, define ^'£{W) = ^>^^(e-"l^ + 
Vl — e~^'^W') for u G [0, +oo[. Using (|4.10p . it is clear that, for 9 < t, we have for u £ [0, +oo), 
< c < ^^^"^{W) < C. Combining the bounds on and ^^^^ yields, for 6' < i and u G [0, +oo), 

< < ^>|^(W)^>^;'(T^) < C^. (4.11) 

Finally, let 

g{z)=l\~^B (B'{{^'z^{W),^^{W))L2iio,T]))\Zt - nZt) = z) du 

= e-"E iyE' (^1^ ^%{W)<!>'^^{W)de^ \Zt - E{Zt) = du. 

The bounds obtained in ()4.1ip immediatly yield < ct < g{z) < Ct. Thus, Proposition 12.21 
concludes the proof of Theorem 14. 6i □ 

Remark 4.7. It is also possible to derive Gaussian density estimates for more complex equations 
than the one dealt with in this section. Hypotheses have to be changed in each case, making it 
difficult to state a general result with reasonable hypotheses covering most cases. 
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